Phase space patterns of quantum transport on ordered and disordered networks 
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In this paper, we consider the quantum-mechanical phase space patterns on ordered and disordered 
networks. For ordered networks in which each node is connected to its 2m nearest neighbors (m on 
either side), the phase space quasi-probability of Wigner function shows various patterns. In the 
long time limit, on even-numbered networks, we find an asymmetric quasi-probability between the 
node and its opposite node. This asymmetry depends on the network parameters and specific phase 
space positions. For disordered networks in which each edge is rewired with probability p > 0, the 
phase space displays regional localization on the initial node. 
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I. INTRODUCTION 

Quantum walk (QW) , which is a generalization of the 
classical random walk, has attracted a great deal of atten- 
tion from the scientific community. The continuous inter- 
est in the study of quantum- mechanical transport process 
can be partly attributed to its broad applications in the 
field of quantum information and computation 
In recent years, two types of quantum walks exist in the 
literature: the discrete-time quantum coined walks and 
continuous- time quantum walks [s', '6*1 . Both discrete- 
time quantum coined walks and continuous-time quan- 
tum walks have been argued to give an algorithmic 
speedup with respect to their classical counterparts [?[. 

In classical physics, the dynamical behavior of a sys- 
tem is described by phase space variables, such as posi- 
tion and momentum. A plot of position and momentum 
variables as a function of time shows the phase diagram 
of classical transport process. In contrast to the clas- 
sical transport, the quantum-mechanical transport hap- 
pens in the Hilbert space P. Such transport process 
needs to be formulated in phase space as a unified pic- 
ture of the classical transport. This can be done by the 
widely used Wigner function, which transforms the wave 
function of a quantum mechanical state into a function 
in the position-momentum space analogously defined in 
the classical phase space [91 . It is shown in Ref. that 
integrating the Wigner function along the lines in phase 
space is a positive value of probability and gives the cor- 
rect marginal distributions. However, the negativity of 
Wigner function provides an indication of non-classical 
behavior. 

The phase space method (Wigner function) provides 
a very useful tool for the study of quantum states in 
the field of statistical mechanics, quantum chemistry, 
molecular dynam ics, scatterin g th eory, quantum optics, 
etc [O, [H, [H [11 [H [H, 111]. There are various 
approaches available to generalize the Wigner function 
for quantum systems with a finite-dimensional space of 
states [18j]. In Ref. [9,], Wootters have introduced a dis- 
crete version of the Wigner function that has all the de- 
sired properties only when TV is a prime number. The 



phase space defined by Wootters is an x grid (A^ 
is prime) and a Cartesian product of such spaces cor- 
responding to prime factors of N in the most general 
case [l^. Notably, a more general discrete Wigner func- 
tion defined for a system with arbitrary values of N 
has been introduced by Hannay and Berry in the stud- 
ies of semiclassical properties of classically chaotic sys- 
tems [2^ . This version of Wigner function is used in sev- 
eral contexts and recently applied to analyze the phase- 
space representation of quantum computers and algo- 
rithms [H, Hi], ll^l- In this case, the phase space is con- 
structed as a grid of 2N x 2N points where the state 
is represented in a redundant manner (only N x N 
them are independent) (loj . Recently, Miilken et al pro- 
pose a version of discrete phase space of Wigner func- 
tion, which is defined for continuous-time quantum walks 
(CTQWs) on a one-dimensional discrete network with 
periodic boundary conditions • This kind of discrete 
Wigner function recovers the correct marginal distribu- 
tions when it is added over the horizonal and vertical 
lines. However, it is not positive when added over the 
general lines [l9| in phase space. This is an unique feature 
differs from the usual Wigner function, which is positive 
when added over any lines . 

Here, we use the version of discret e p hase space of 
Wigner function proposed by Miilken In Ref. p^ . 
Miilken et al formulate CTQWs in phase space on 
a network of size N whose nodes are enumerated as 
0, 1, A^ — 1. The Wigner function has the form of a 
Fourier transform as follows llOll. 



W(x,k,t) 



1 

TV ^ 



e « <x-y\p{t)\x + y>, (1) 



Where p{t) is the density operator of a pure state, k and 
X denote the phase space coordinate of positions. The 
summation over y in the interval [0, A^ — 1) can be car- 
ried out in any A^ consecutive values. Considering an 
initial exciton begins at node j, the time evolution of the 
associated state \ j > is given by |j, t >= e~'^'|j >. Sup- 
pose En and |g„ > are the nth eigenvalue and eigenstate 
of the Hamiltonian of Laplace matrix, the time indepen- 
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dent Schrodinger equation is H\qn >= En\qn >, where 
\qn > spans the whole accessible Hilbert space and forms 
an orthonormal complete basis set, i.e., < qn\qi >= Sni, 
X^n 1 9" 9"! ~ ^- Inserting the complete set condition 
of the eigenstates into the time evolution equation, we 
get, 

\j,t>^ e-^"'\j ^e-*^"*|g„ X qn\j > . (2) 



The density operator of the system is p{t) = \j, t >< 
t, j\. Substituting the above Equation into the Wigner 
function, we have, 



• < a; - y\qn >< qn\j >< j\qi >< qi\x + y > . 



In this paper, we use the above Equation to consider 
phase space patterns of CTQWs on ordered and disor- 
dered networks. For ordered networks, the topology or- 
ganizes in a very regular manner, i.e., each node of the 
network is connected to its 2m nearest neighbors (m on 
either side). For disordered networks, we employ the 
famous WS model [1^, which triggers a surge of re- 
search of small-world networks in the field of complex 
networks [H, [H [H [H . In the WS model [H, each 
connection of the regular networks is rewired with prob- 
ability p. Tuning the rewiring probability p interpolates 
the network topology between order (p = 0) and disor- 
der (or random with p — 1). The intermediate value 
< p < 1 corresponds to the small- world region [2^ . 

The paper is organized as follows. In Sec. II, we con- 
sider the phase space patterns on ordered networks which 
correspond to the WS model with rewiring probability 
p = 0. In Sec. Ill, we consider the phase space patterns 
on disordered networks. Conclusions and discussions are 
given in Sec. IV. 



A. Bloch ansatz and Wigner function 

The Hamiltonian (H) of the system for CTQWs is re- 
lated to the Laplace matrix (A) of the connected net- 
works a.s H = J A. Here, for the sake of simplicity, we 
assume the transmission rate 7 for all the connections to 
be equal. The nondiagonal elements Aij equal to —1 if 
nodes i and j are connected and otherwise. The di- 
agonal elements An equal to the number of total links 
connected to node i, i.e.. An equals to the degree of node 
i. Therefore, the Laplace matrix A of ordered networks 
takes the following form. 




J, 



if i 

if i — j ± z £ 
Otherwise. 



[l,m] 



(4) 



The Hamiltonian acting on the state \ j > can be written 
as 

m 

H\j >= (2m -I- l)|j > — \j + z >, z £ Integers. 

z——rn 

(5) 

The above Equation is the discrete version of the Hamil- 
tonian for a free particle moving on the network. Using 
the approach of Bloch function [3J| in solid state physics, 
the time independent Schrodinger equation reads 



H\lPn >= En\lpn > 



(6) 



The Bloch states \4>n > can be expanded as a linear com- 
bination of the states \j > localized at node j [s^l, 



N-l 



\j> 



Substituting Eqs. © and ^ into Eq. 
eigenvalues (or energy) of the system. 



(7) 



we obtain the 



II. PHASE SPACE PATTERNS ON ORDERED 
NETWORKS 



In this section, we consider the quantum-mechanical 
phase space on one-dimension regular networks of N 
nodes in which each node is connected to its 2m near- 
est neighbors (m on either side). This generalized reg- 
ular network has broad applications in various coupled 
dynamical systems, including biological oscillators [29|, 
Josephson junction arrays [30], neural networks 31 1, syn- 
chronization [3^ . small- world networks [33] and many 
other self-organizing systems. We compute the phase 
space distribution on such general network with peri- 
odic boundary conditions in the framework of Bloch 
ansatz [3^ . which is commonly used in solid state 
physics. 



En ~ 2m — 2^ cos(j0„). 



(8) 



The periodic boundary condition for the network requires 
that the projection of IV'n > on the state |iV > equals to 
that on the state |0 >, thus 0„ — 2nTT/N with n integer 
and n S [0, iV). Replacing |g„ > by the Bloch states 
\ipn > in Eq. ([3]), we can get the Wigner function as 
follows, 



Wj{x, k,t)^^ Ef=o^ " Y.n,i e" 

.gie„ (x-y-j) ^ie, (j-x-y) 



it{E„-Ei) 



2^y=o e 



(9) 

The summation over y can be written as Af(5(jv+fc-ri-i)j 
where S^j^^^-n-i) takes value 1 if {N + k — n — l) equals to 
(or mod N) and otherwise. Thus the Wigner function 



FIG. 1: (Color online) Wigner functions Wj{x,k,t) on net- 
works of size A'' = 50 with m — 1 ((row 1), m = 2 (row 2) and 
m = 3 (row 3) at times t — 1,2, 5, 20 (columns (a)-(d)). The 
initial node is at j = N/2 — 25. Red regions denote positive 
values of Wj{x,k,t), blue regions denote negative values and 
white regions close to value 0. 

can be simplified as, 

71=0 

(10) 

Substituting the expression of eigenvalues of Eq. ^ into 
the above equation, we obtain discrete Wigner function 
of CTQWs on one-dimension ordered networks. 

B. Time evolution of the Wigner function 

We consider the time evolution of quantum-mechanical 
phase space according to Eq. ^TU\i . Figure [1] shows a con- 
tour plot of the Wigner function on ordered networks of 
size iV = 50 with different values of m at different times. 
At very small time scales, the Wigner function is localized 
on the stripe at the initial node and its opposite node. 
As time increases, the Wigner function spreads over the 
whole network. On short time scales, the Wigner func- 
tion has a very regular structure until the wave fronts 
between the initial node and opposite node start to in- 
terferes with each other. Interestingly, the phase space 
structure is more complex on highly connected networks 
compared to that on the cycle network (to =1), and the 
nodes are populated more quickly on highly connected 
networks. 

In Fig. [2 we plot the Wi gner function of CTQWs on 
networks of size = 51. The phase space structure is 
quite similar to that on even- numbered networks. We 
note that the phase space pictures provide us more infor- 
mation of the underlying dynamics than the transition 
probabilities. The phase space patterns alternate with 
time frequently. At long time scales, the phase space be- 



FIG. 2: (Color online) Corresponding plots of Fig. [l]for A'^ = 
51. 



comes irregular and we find that the structure of phase 
space has more regularities on even-numbered networks 
than that on odd-numbered networks with the same value 
of TO and t. This indicates the higher topological sym- 
metry of even-numbered networks. 

It is worth noting that the pattern of Wigner function 
is symmetrical about the axis for to = 1. As 771 increases, 
such behavior disappears and Wj{x, k, t) displays central 
symmetry at the phase space center x = k = j . At the 
initial time i = 0, on even- numbered networks (or odd- 
numbered networks) , the patterns of Wigner function are 
the same for different values of to. This can be concluded 
from Eq. pO)) . For the case of even TV, Wj{j, k, 0) equals 
to 1/A^ for arbitrary k. At the opposite node x — j+N/2, 
Wj {j + N/2, k, 0) equals to 1 /N for even fc and - 1 /N for 
odd k. The nonzero values of the Wigner function at 
X — j + N/2 continue to show up at later times. The 
analysis for odd N is similar but the patterns are differ- 
ent. The Wigner function Wj{x,k,0) equals to 1/A^ for 
X = j and otherwise. The nonzero of Wj{x, k,0) at 
the opposite node x = j + N/2 {N € Evens) is a natural 
consequence of the periodic boundary conditions of the 
regular networks [l3|- At later times, Wj{x, k, 0) involves 
contribution from all the eigenstates and the nodes con- 
nected to the excitation node get populated, resulting the 
semicircle-like areas on highly connected networks. We 
remark that during all the time the patterns show central 
symmetry at x = fc = j, to this end, we conjecture that 
central symmetry is an intrinsical feature of the phase 
space patterns. 



C. Long time averages 

The Wigner function of a specific position {x, k) fiuc- 
tuates around a constant value, thus it is interesting to 
study the long time averaged phase space patterns. The 
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time limiting Wigner function is defined as, 

1 

>Vj(x,fc)= lim - / Wj{x,k,t)dt. (11) 

For the ordered networks, the hmiting Wigner function 
can be simphfied as. 



N~l 
n=0 



(12) 



This expression shortens the numerical time of compu- 
tation considerably compared to Eq. (jll[) . Therefore, we 
consider the long time averaged phase space structure 
according to Eq. (fT2|) . 

For the cycle network (m = 1) in which each node is 
only connected to its two nearest neighbors, the limiting 
phase space has a simple structure. For even-numbered 
networks {N e Evens), the limiting phase space struc- 
ture is 

k ^ Q,k E Evens and arbitrary x, 
WUx,k)^{ 1/N, k = andx^j,j±N/2, 



0, 



eLsewnere. 



(13) 

If the network size N is an odd number, we can also ob- 
tain the limiting Wigner function according to Eq. (fT2l). 
which is summarized as, 

{1/iV^, k ^ and arbitrary x, 
l/N, k^Oandx^j, (14) 
0, elsewhere. 

which confirms the results in Ref. [s^. 

For other values of m, the limiting phase space distri- 
butions can also be calculated in the same way, but such 
process is complicated for large value of m because of 
the nontrivial degeneracy distribution of the eigenvalues. 
Here, we report the numerical results of phase space pat- 
terns on highly connected networks according to Eq. (H^). 

Fig. [3] shows the phase space patterns on networks of 
= 50 with different values of m. We note that the 
phase space structure of m = 1 is the same as the phase 
space structure of m = 3. After a careful examination, 
we find that m = 1 and m = 3 have the same degeneracy 
distributions of the eigenvalues. This explains the same 
phase space patterns of m = 1 and ni — 3. For m = 2 
and m = 4, there are some significant stripes in the phase 
space. Such nontrivial stripes reflect the topological sym- 
metry of the considered networks. On the contrary, the 
phase space patterns on networks of A'^ = 51 with to = 2, 
TO = 3 and to = 4 are the same as the structure with 
TO = 1, which is described in Eq. (jl4p . One may conjec- 
ture that the phase space patterns do not alter when in- 
creasing the connectivity to on odd-numbered networks. 
But this is not true for some particular value of network 
size N. For instance, on a network of A'^ = 75 and to = 2, 




FIG. 3: (Coloronline) Long time limiting Wigner function 
Wj{x,k) on ordered networks of size = 50 for m = 1 (a), 
m = 2 (b), m = 3 (c) and m — 4 (d). The initial node 
is at j — N/2 — 25. Red blocks denote positive values of 
yVj{x,k), blue blocks negative values, and the white regions 
denote zero. 



we find significant stripes at some phase space positions 
(See Fig. HI). 

The patterns of Wj {x, k) in Fig. [3] are the same for odd 
k and some even k. However, for some particular values 
of TO and even fc, the patterns are quite different (See the 
stripes at fc = 10, 20, 30 and 40 in Fig.[3i;b) and (d)). Ac- 
cording to Eq. (|12p . Wj{x, k) depends on the degeneracy 
of the eigenstates, i.e., the Kronecker symbol 5E„,EN+k-n- 
For fc = 0, ^£;„,£;jv+fe-,i equals to 1 for all the values of n, 
thus the sum in Eq. (|12p contains A^ exponential terms, 
which results in yVj(a;,0) = 1/A^ for x — j,j ± N/2 and 
Wj{x,0) — otherwise. When k is odd, SE„,EN+k-,t — 
holds for all the values of n. For some values of even fc. 



^ equals to 1 when n = k/2 and n = N/2+k/2, 
this leads to yVj(a;, /c) = 2/A^^ {k e arbitrary x). How- 
ever, for some other values of even fc, SEr^.EN+k-n does not 
vanish and Wj{x, k) is oscillatory on the horizontal lines 
located at A: = 10,20,30 and 40. The oscillating strip 
is caused by the interference between the positive strip 
and the mirror image. This is similar to the case of deco- 
herence of quantum walks in phase space (isj . in which 
there are interference fringes in lines between the two po- 
sition eigenstates, and all the vertical lines have their cor- 
responding oscillatory counterparts originated from the 
boundary conditions [isj . 

In Fig. [3l we find a symmetric structure of the phase 
space distributions in both the x and k directions. That 
is to say, Wjix, k) equals to Wj{x+N/2, k) {x G [0, N/2)) 
for all the values of k; Wj(x, k) equals to Wj(a;, k + N/2) 
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FIG. 4: (Coloronline) Long time averaged phase space pat- 
terns on a network of A'^ = 75 and m = 2. The initial node 
is at j = 37. Red blocks denote positive values of yVj{x, k), 
blue blocks negative values, and the white regions denote the 
value of 0. 
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FIG. 5: (a)Contour plot of Wj{x,k + N/2) - yVj{x,k) on 
networks of size A'' = 60 with m = 2 (row 1), m = 3 (row 
2) and m = 4 (row 3). (b): the same as column (a) but 
for Wjix + N/2,k) - yVj{x,k) (x € [0,iV/2)). Red blocks 
denote positive values of the discrepancy, blue blocks negative 
values, and the white regions denote value of (symmetric 
quasi-probability in the phase space). The last two columns 
show marginal distributions of Wj{x, k) (column (c)) and 
X]fc Wj(a;, k) (column (d)) for m = 2, m = 3 and m = 4 (rows 
1-3). The initial node is at j = N/2 = 30. 



[k €E [0, iV/2)) for all the values of x. Such symmetric 
phase space structure exists on some even-numbered net- 
works. However, for some even- numbered networks with 
certain values of to, this is not true. For instance, on 
networks of size = 60 with m = 2, to = 3 and m = 4, 
the Wigner function Wjlx, k) differs from Wj{x + N/2, k) 
{x € [0, N/2)) for some values of k and Wj{x, k) also dif- 
fers from Wj{x, k+N/2) [k g [0, iV/2)) for some values of 
X. Such asymmetry depends on the specific phase space 
positions {x,k). In Ref. [s^, the authors find an asym- 
metry of transition probabilities for the starting node and 
its mirror node on two-dimensional networks, their defini- 
tion of mirror node is based on geometry symmetry of the 
network. In this paper, we define the mirror node a; of a 
given node x to be its opposite node, i.e., x = x + N/2 
(or k = k + N/2 in the longitudinal direction). We study 
symmetric and asymmetric structure of the phase space 
on the X and k direction separately. The asymmetric 
phase space is particularly characterized by the differ- 
ence between Wjix, k) and Wjix + N/2, k) {x G [0, N/2)) 
in the x direction, and difference between yVj{x,k) and 
yVj{x,k + N/2) (k G [0,N/2)) in the k direction. There- 
fore, we use the quantities Wjix + N/2,k) — Wj{x,k) 
and yVj (x, A: + N/2) — >V, (a;, /c) to detect the asymme- 
try of the phase space. Columns (a) and (b) of Fig. O 
show these two quantities in the phase space. The non- 
white blocks in the plot indicate asymmetric phase space 
positions. From the figure we find that the asymmetric 
structure in the phase space is very complex. The phase 
space positions in which asymmetry occurs is dependant 
on the precise values of network parameters N and m, as 
well as the specific phase space positions. 

It is interesting to note that, in Figs. [T] and O Wj{x, k) 
has more identical values than Wj{x,k,t) in the phase 
space plane. Wj{x,k,t) and yVj{x,k) display central 
symmetry and mirror symmetry respectively. This sug- 
gests that the limiting Wigner function >Vj(x,fc) has 
higher symmetry in the structure compared to the in- 
stantaneous Wigner function Wj{x,k,t). 

D. Marginal distributions 

Summing the Wigner function along lines in phase 
space gives the marginal distributions. Of course, such a 
process loses the detailed information of the phase space. 
Here, we consider marginal distributions of the long time 
averaged Wigner function summing over the x direction 
and k direction respectively. Summing over all k recovers 
the time limiting transition probabilities, i.e., 

k 

Columns (c) and (d) of Fig. [5] show the marginal distri- 
butions of the sum J^x'^ji^^^) ^^"^ J^k'^ji^^^)- 
find that the marginal distributions obtained by summing 
over X for m = 2, to = 3 and to = 4 are the same (See 
column (c) in Fig. [5]). equals to 2/N for 
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FIG. 6; (Color online) Quantity x^+Jv/aj ~ X^,i a func- 
tion of X on networks of = 60 with different values of m. 
The initial node is at j = N/2 = 30. The nonzero value of 
{Xx+N/2,j — Xx,j) presents asymmetric probabilities between 
Xx+N/2,j and Xx,j- 




FIG. 7: (Color online)Phase space patterns on networks of 
= 50 with exponential disorder. The three rows are for 
m — 1 (row 1), m = 2 (row 2) and m = 3 (row 3) while the 
four columns correspond to the A — 0.5 (column (a)), A — 0.6 
(column (b)), A = 0.7 (column (c)) and A = 0.8 (column (d)). 
The initial node is at j = N/2 = 25. Red regions denote 
positive values of Wj{x, k), blue regions negative values and 
white regions denote values close to 0. 



even k and for odd k. In contrast, the marginal distri- 
butions Xx.j for different values of m are quite different 
(See column (d) in Fig. [5]). Interestingly, there is a large 
probability to be still or again at the initial node (a; — j) 
and at the opposite node {x = j — N/2). In analogous to 
the analysis of the phase space asymmetry, we consider 
asymmetry of the limiting transition probability. We find 
that Xx+N/2.j differs from Xx,j{x G [0, 7V/2)) for some 
particular values of x on networks with certain values of 
N and m. Fig. [H] shows the quantity Xx+N/2,j — Xx.j for 
all the values of x in the range [0, A^/2). The nonzero 
values Xx+N/2.j ~ Xx,j in the plot indicate asymmetric 
transition probabilities at the corresponding nodes. 

As shown in Ref. [36j , the asymmetries of transi- 
tion probabilities originate from different contribution of 
eigenvalues to Xx,j- It turns out that there arc more 
contribution to X x.j in the asymmetric cases than in the 
symmetric cases 36|. The argumentation for our case is 
similar. Although all the eigenvalues can be analytically 
obtained, a complete analysis of all possible differences of 
eigenvalues requires extensive work and is clearly beyond 
the scope of this paper. 



connections, thus does not change the topology of the 
ordered networks. On the contrary for networks with 
topological disorder in which each connection is rewired 
with probability p > 0, the topology of the networks be- 
comes irregular. 

In this paper, we focus on the influence of disorder on 
the long time averaged phase space structure. Since the 
analytical expressions for ordered networks do not apply 
any more, we numerically calculate the Wigner function 
of CTQWs using the software Mathematica. In order to 
save the computational time of the numerical integrals, 
we combine Eqs. ([3]) and (fTTjl as follows, 

■ <x~ y\qn X qn\j >< j\qi >< qi\x + y > . 

In the following, we use the above Equation to report the 
numerical results of Wigner function on networks with 
disorder. 



A. static disorder 



III. PHASE SPACE PATTERNS ON 
DISORDERED NETWORKS 

In this section, we study the phase space patterns on 
networks in the presence of two kinds of disorder: static 
disorder and topological disorder. We import static dis- 
order by adding a perturbed Hamiltonian A to the origi- 
nal Hamiltonian H as done in Ref. [37[ . This assumption 
do not create new connections or remove the existing 



Analogous to the method in Ref. [37|, we consider 
the exponential diagonal disorder whose off-diagonal el- 
ements of the perturbed Hamilton A are and diagonal 
elements follow an exponential function as, 

A„- = e^^, (17) 

where A is the disorder parameter. The exponential dis- 
order in the above Equation has the advantage that the 
perturbed Hamilton A possess a simple form in the Bloch 
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FIG. 8: (Color online)The same plots as Fig.[7]but for N = 51. 



representation, and analytical solutions may be possible 
using the perturbation theory. Here, we only give numer- 
ical results for the phase space structure for such disorder. 
Figs. [7| and [8] show the phase space structure on networks 
of size = 50 and = 51 for different values of m and A. 
The first column of the figures shows Wj{x, k) for A = 0.5 
on networks with m = 1, to = 2 and m = 3 (rows 1-3). 
For this weak disorder, the phase space patterns have a 
very strange structure and the patterns of the ordered 
networks are destroyed. Increasing the disorder parame- 
ter A, the patterns change drastically. The phase space 
structure gets suppressed and a localized region forms at 
the initial node x = j. In contrast to the patterns in 
Ref. 37], the Wigner function here has negative values 
on highly connected networks with large disorder (Com- 
pare the plots in the last column to plots in Ref. [STj). 
The patterns in phase space on odd-numbered networks 
(Fig. [8]) have a similar structure as the even-numbered 
networks. However, differences are also visible and it 
seems that there are more stripes on even-numbered net- 
works (Compare the corresponding plots in Fig. [7] and 
Fig.©. 

Compared to the patterns in Ref. [37] , the phase space 
structure for networks with exponential disorder shows a 
strange patterns. Such difference is induced by the dis- 
tinct type of disorder. The perturbed Hamiltonian in 
Ref. f37| has a Gauss form, such disorder leads to a sym- 
metrical phase space pattern in the x—k plane. Here, the 
exponential disorder has a heterogeneous strength, i.e., 
sites labeled as large numbers have large exponential dis- 
order (See Eq. P?|) ). this heterogeneous disorder results 
in the complex pattern of the observed phase space. For 
larger disorder, the difference of phase space between net- 
works with differen connectivity becomes invisible. We 
believe for sufhcient strong disorder, the phase space dis- 
plays a complete localization at the initial node. 



FIG. 9: (Color online) Phase space patterns on WS networks 
of A'^ = 50 with m = 1 (row 1), m — 2 (row 2) and m = 
3 (row 3). The four columns correspond to different values 
of disorder parameter p: p — 0.05 (column (a)), p = 0.1 
(column (b)), p — 0.2 (column (c)) and p = 0.5 (column (d)). 
The initial node is at j = N/2 = 25 and all the plots are 
averaged over 200 realizations. Red regions denote positive 
values of Wj{x,k), blue regions negative values and white 
regions denote values close to 0. The colormaps are the same 
for all the plots. 



B. topological disorder 

Since the topological structure of the WS model is not 
single, we average the time limiting Wigner function over 
distinct realizations. The ensemble average of the Wigner 
function provides us a holistic view on the phase space 
structure of disordered networks. Figs. [9] and [10] show the 
ensemble averaged phase space patterns on disorder net- 
works of size = 50 and A^ = 51 with different values of 
TO and p. It is found that the stripes on regular networks 
are destroyed and localized regions form around the ini- 
tial node X = j (Compare Fig. [3] and Fig. [5]). Increasing 
the rewiring probability p (disorder parameter in the WS 
model) , the patterns change profoundly. The phase space 
structure gathers together at the initial node x = j for all 
k on networks with large disorder. However, the stripes 
on ordered networks are still visible even for large disor- 
der. When the disorder parameter p is fixed, the effect of 
localization becomes weak on networks with more con- 
nectivity (Compare the plots in the same columns). For 
odd-numbered networks, the patterns of localized central 
regions are rather comparable to that on even-numbered 
networks. The other regions in the phase space are anal- 
ogous to the patterns of the corresponding ordered net- 
works, although the local patterns around the a; « j are 
changed. 

Summing the averaged limiting Wigner function over 
k gives the averaged transition probability < Xx,j >= 
J2k < Wj{x,k) >. Fig. [TT] shows this marginal distri- 
bution < Xx,j > for even and odd A^ with different con- 
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FIG. 10: (Color online)The same plots as Fig. [9] but for N ■ 
51. 



nectivity and disorder strength. In the figure, we can see 
that there is a remarkable localization at the initial node 
when m — 1. As m increases, the localizations become 
weak and nearly the same for m = 2 and m = 3 (Com- 
pare the corresponding dots and curves in row 2 and 3). 
Interestingly, we note that the marginal distributions on 
even and odd numbered networks are nearly the same 
and there is only one peak a.t x ~ j. This feature differs 
from the case on ordered networks where there are two 
peaks {x = j and x = j + N/2) for even N and only one 
peak (x = j) for odd N. The two peak marginal distribu- 
tion on regular even-numbered networks is a consequence 
of the symmetry of the network topology. When disor- 
der was added, such topological symmetry is destroyed, 
resulting in the single-peaked marginal distributions on 
disordered networks. 



IV. CONCLUSIONS AND DISCUSSIONS 

In conclusions, we have studied the quantum- 
mechanical phase space patterns on ordered and disor- 
dered networks. On ordered networks where each node 
connects to its 2m nearest neighbors (m on either side), 
the phase space quasi-probability of Wigner function 
shows various patterns. In the long time limit, the phase 
space structure presents different kinds of stripes. Such 
striped patterns are related to the specific network size N 
and connectivity parameter m. Interestingly, if the net- 
work size N is an even number, we find an asymmetric 
quasi-probability and transition probability between the 
node and its opposite node. This asymmetry depends 
on the network parameters and specific phase space po- 
sitions. On disordered networks in which each edge is 
rewired with probability p > 0, the phase space displays 
regional localization on the initial node. 

The asymmetry of the quasi-probability and transition 
probability is a novel phenomenon, which does not exist 
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FIG. 11: (Color online) Transition probability < Xx,j >= 
"^f. < yVj{x,k) > on networks of size A' = 50 (column (a)) 
and A = 51 (column (b)) with m = 1, m = 2 and m = 3 (rows 
1-3). In each figure, the rewiring probability p takes values 
0.05 (squares), 0.1 (circles), 0.2 (triangles) and 0.5 (rhombus). 



in the cycle graph with m = 1. However, we are un- 
able to predict which particular parameters of N and m 
or which phase space positions (x,j) are related to such 
asymmetry. Is there relation between the phase space 
asymmetry and transition probability asymmetry? Such 
question is interesting and requires a further study. The 
phase space patterns on disordered networks suggest that 
there are localizations in phase space on small- world net- 
works. Although both the static disorder and topologi- 
cal disorder lead to localizations at the initial node, lo- 
calizations on small-world networks indicates that it is a 
generic property of disordered systems and has important 
consequences for quantum walk algorithms and quantum 
communication 1381. l39l . 
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